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ABSTRACT :  The  paper  first  presents  a  review  of  prior  developed  methods  for  the 
determination  of  deflections  of  the  vertical  from  astrogeodetlc  and  Inertlally 
obtained  data  under  consideration  of  gyro  bias  eliminations  In  Lltton's  local- 
level  system.  It  then  derives  a  Wiener-type  optimal  solution  In  semi-flat  ter¬ 
rain  under  utilization  of  the  collocation  method  In  physical  geodesy  with  prior 
emphasis  on  non-linear  gyro  bias  removal  and  discusses  the  ramifications  of  pre¬ 
sently  and  potentially  available  gyro  and  accelerometer  hardware.  Thereafter,  a 
refined  advanced  method  and  deflection  determination  In  strongly  mountainous  ter¬ 
rain  Is  addressed  under  possible  modification  of  the  semi-flat  terrain  solutions 
and  employment  of  spatial  collocation.  Finally,  consideration  Is  given  to  the 
associated  problem  of  area  adjustment  solutions. 

1.  INTRODUCTION.  The  selective  Interpolation  of  deflections  of  the  vertical  by 
means  of  Initial  and  terminal  astrogeodetlc,  and  inertial  data  has  been  pursued 
by  the  U.S.  Army  Engineer  Topographic  Laboratories  (ETL)  and  the  Geodetic  Survey 
of  Canada  since  about  1976.  The  Inertial  equipment  employed  has  been  the  Rapid 
Geodetic  Survey  System  (RGSS),  developed  by  Litton  Systems  for  ETL.  Presently 
installed  gyros,  accelerometers,  and  velocity  quantizers  have  permitted  average 
deflection  component  accuracies  of  1.5  arcsec  rms  for  single  runs  of  60  km  length 
or  approximately  2  hours.'  The  Installment  of  available  superior  hardware  In  the 
horizontal  channels  of  the  RGSS  and  the  utilization  of  Improved  post-mission  data 
reduction  methods  Is  expected  to  result  In  average  deflection  component  accuracies 
of  about  0.4  arcsec  rms  for  single  runs.  Network  adjustments  under  consideration 
of  data  obtained  from  quasi-parallel  runs  may  thus  facilitate  average  deflection 
component  accuracies  of  0.3  arcsec  rms.  In  this  respect.  It  should  be  emphasized 
that  the  combination  of  new  gyros  and  accelerometers  with  relatively  small  correla¬ 
tion  times  and  Improved  data  reduction  methods  Is  particularly  effective.  This  pa¬ 
per  reviews  In  the  next  section  previously  developed  methods  for  the  determination 
of  deflections  of  the  vertical.  In  the  third  section,  a  Wiener-type  solution  for 
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semi -flat  terrain  Is  presented  together  with  ramifications  of  present  and  Im¬ 
proved  gyro  and  accelerometer  hardware.  The  fourth  section  outlines  a  refine¬ 
ment  of  the  optimized  method.  The  Influence  of  mountainous  terrain  Is  consid¬ 
ered  In  the  fifth  section,  and  the  sixth  section  addresses  the  problem  of  ad¬ 
justment  solutions. 

2.  REVIEW  OF  PREVIOUSLY  DEVELOPED  METHODS:  Baussus  von  Luetzow  [1981a]  pre¬ 
sented  In  detail  the  approximate  mathematical  solution  for  deflections  of  the 
vertical  as  presently  applied  In  the  context  of  RGSS  utilization.  In  this  res¬ 
pect,  the  basic  solution  for  the  prime  deflection  of  the  vertical  Is 

'i ■  <«. ♦  V k >  - ii 

18  8  tn  g 

t ‘I  -tn 

+  /  (y  -  Y>dt  J  (Y-?)dt  (1) 

O  tn  O 

In  eq.  (1),  £  Is  the  prime  deflection,  ¥E  Is  the  systematic  platform  tilt  about 
the  east  axis,  x  Is  the  static  north  accelerometer  measurement,  aN  Is  the  cor¬ 
related  east  accelerometer  error,  y  is  the  east  axis  angular  drift  rate  error, 

Y  Is  the  associated  constant  gyro  bias,  g  Is  normal  gravity,  and  t  Is  time.  The 
subscripts  refer  to  Initial  time  t0,  time  at  the  start  of  a  vehicle  stop  with 
accelerometer  reading  tA,  and  terminal  mission  time  tn: 

Equation  (1)  may  be  supplemented  by 

«<1  -“«>♦£  (2) 

to  account  for  astrogeodetlc  deflection  errors  and  accelerometer  bias  errors. 

For  a  straight  traverse,  the  last  term  In  eq.  (2)  tends  to  cancel  out. 


The  determination  of  *E  and  associated  gyro  bias  functions  Is  accomplished 
by  means  of  the  equations 

**  ■  1  (if  ♦  «  y*  -«  (3) 

1 4  j  \  w 

♦  •  )  (5) 

where  ?N  is  the  systematic  platform  tilt  about  the  north  axis,  ?  Is  the  system¬ 
atic  azimuth  platform  attitude,  a  is  the  constant  azimuth  axis  angular  drift  rate 
bias,  e  is  the  constant  north  angular  drift  rate  bias,  nN  is  the  north  earth  rate, 
and  J)  is  the  vertical  earth  rate. 

I* 

■  C  estimated  ^  "  *  observed  (6) 

■  n  estimated  *  11  observed  (7) 

m  *  estimated  *  *  observed  (8) 

It  should  be  emphasized  that  the  errors  relating  to  the  computation  of  and 
the  static  accelerometer  reading  In  eq.  (1)  are  reflected  In  the  last  four  terms. 
Examination  of  equations  (3)  -  (5)  r eve} Is  that  In  moderate  latitudes  there  exists 
relatively  small  coupling  for  time  Intervals  not  exceeding  2  hours.  For  this  rea¬ 
son,  linear  approximations  »E  »  yt  and  FN®»(Ft  have  been  used  with  success.  High¬ 
ly  accurate  deflection  determinations,  accomplished  by  means  of  an  advanced  RGSS, 
would,  however,  require  consideration  of  a  terminal  azimuth  error,  modified  by  a 
periodically  applied  Kalman  filter  correction. 

The  rms  deflection  error  o^[t±)  can  be  computed  by  covariance  analysis  Involving 
the  terms  without  parentheses  In  eq.  (1).  Under  Inclusion  of  the  first  two  terms 
of  eq.  (2)  It  Is 

var  •  var  a*  ♦  var  Trt  ♦  (1  -  ||)*  var  t0  ♦  v*r  ^  (9) 
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where  var  a1  Is  the  accelerometer-induced  variance  and  var  designates  the 
gyro-induced  variance. 

The  use  of  static  accelerometer  measurements  In  conjunction  with  Initial 
and  terminal  deflections  together  with  a  simplified  Kalman  filter  originated 
by  Huddle  [1977].  A  weighted  least-squares  solution  for  deflections  of  the 
vertical  was  developed  by  Baussus  von  Luetzow  [1978]  under  Integration  of  the 
pertinent  system  of  differential  equations,  constant  travel  Intervals  between 
vehicle  stops,  utilization  of  a  limited  number  of  deflection  components  together 
with  collocation  estimation,  and  velocity  observations.  In  this  approach,  gyro 
biases  were  treated  as  random  variables.  For  practical  purposes.  It  Is  necessary 
to  record  the  total  velocity  history  or  to  drive  the  RGSS  at  a  constant  velocity. 

3.  ADVANCED  METHOD  FOR  SEMI-FLAT  TERRAIN.  Since  the  determination  of  ac¬ 
cording  to  eq.  (1)  Is  Impaired  because  of  the  presence  of  four  noise  terms,  a 
Wiener-type  solution  under  consideration  of  all  or  of  adjacent  Cj-data  in  con¬ 
junction  with  a  signal  covariance  function  can  be  formulated  as  shown  by  Baussus 
von  Luetzow  [1981a].  For  this  purpose,  eq.  (1)  under  Inclusion  of  the  first  two 
terms  of  eq.  (2)  may  be  formulated  as 

E*  ■  +  n±  (10) 

where  -  E0  +  *Ei  -  -21  Is  a  message  variable,  £±  Is  a  signal  variable,  and  the 
remaining  terms  denote  a  noise  variable  -  n±.  The  collocation  method  In  physical 
geodesy  In  semi-flat  terrain  permits  the  estimation. 

^  >l^i4"1)  ■  *i(Si  4  *i>  (11) 

where  Aa  Is  the  matrix  of  regression  coefficients  at  to  be  computed  and  Is  the 
corresponding  E±  -  matrix.  It  Is  then  In  matrix  form,  with  bars  Indicating  covarl- 


The  solution  for  the  regression  coefficient  matrix  follows  as 

*  "l V  1  (13) 

The  pertinent  noise  parameters  applicable  to  the  present  RGSS  are  a  stand¬ 
ard  deviation  of  0.002*  hr-1  and  a  correlation  time  of  3  hours  for  the  G-300 
gyros  and  r  standard  deviation  of  10  mgal  and  a  correlation  time  of  40  minutes 
for  the  A-200  accelerometers.  Due  to  the  sizable  correlation  times,  the  ad¬ 
vanced  method  does  not  provide  significantly  better  deflection  estimates  than 
the  basic  method.  Considerably  Improved  estimates  would,  however,  result  under 
utilization  of  Intermediate  deflection  constraints.  An  advanced  RGSS,  with  a 
standard  deviation  of  0.0002°  hr-1  and  a  correlation  time  of  about  5  minutes  for 
G300-G2  gyros  and  a  standard  deviation  of  1  mgal  and  an  approximate  correlation 
time  of  5  minutes,  would  generate  data  commensurate  with  the  potential  of  the 
advanced  estimation  method.  In  addition.  It  constitutes  a  statistical  framework 
for  an  area  adjustment  under  utilization  of  data  relating  to  several  traverses. 

4.  REFINEMENT  OF  ADVANCED  METHOD.  The  basic  method  for  the  determination  of 
deflections  of  the  vertical,  presented  In  section  2,  does  not  consider  the  Inter¬ 
action  of  the  two  horizontal  error  differential  equations  of  motion  and  the  three 
gyro  error  differential  equations.  The  latter  are  employed  for  the  determination 
of  constant  gyro  biases.  For  short  traverse  times  and  In  connection  with  present 
RGSS  data  the  basic  method  Is  certainly  adequate.  More  accurate  deflection  deter¬ 
mination  by  means  of  an  advanced  RGSS  may  require  a  refinement  of  the  advanced 
method.  This  Is  consistent  with  observations  made  by  Schwarz  [1980]. 

A  full  refinement  of  the  advanced  method  requires  the  Integration  of  the 
whole  system  of  differential  equations  and  the  application  of  RGSS  Kalman  filter 
corrections  at  stops,  preferably  also  extended  to  #  .  For  practical  purposes, 
vehicle  stop  Intervals  should  be  3  minutes  and  the  vehicle  speed  should  be 
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approximately  constant,  preferably  about  30  km  hr-1.  As  a  result,  5  In  eq.  (1) 
would  be  replaced  by 

5ri  =  +  Ci  +  04) 

V  V 

where  C±  denotes  a  systematic  correction,  alv  and  blv  are  weight  factors  obtained 
by  numerical  integration  under  consideration  of  Kalman  filter  corrections,  and 
n1(2)  is  a  composite  representation  of  modified  gyro  and  accelerometer  errors. 
Further,  £v  and  nv  indicate  ?(ty)  and  ntt^),  respectively.  Equation  (14)  may  be 
formulated  as 

(£aiv^v+  S^iv^v)  a  (50  +  *£.  -  -  n^(2)  (15) 

v  v  1  a 

The  message  estimator  is  then 

V  •• 

Li  =  L±  +  n^2)  =  (50  +  $Ei  -  +  C±)  (16) 

In  analogy  with  equation  (11),  any  signal  variable  may  be  estimated  according  to 

le  =  bi(Li  +  ni(2))  -  h  i  [4]  +  N±(2)|  (17) 

where  is  the  matrix  of  regression  coefficients  and  [l±]  and  N^2)  are  corres¬ 
ponding  signal  and  noise  matrices.  In  accordance  with  eq.  (13),  the  solution  for 
the  regression  coefficient  matrix  follows  as 

Bi  =  led k)  {  (L±J  [£fcl  +  N±(2)Nk(2)}‘  (18) 

In  general,  eq.  (18)  requires  the  coordinates  x(tv),  y(tv),  x±,  and  y±  as  computer 
program  inputs. 

It  is  apparent  from  the  above  analysis  that  the  basic  deflection  estimation 
method  is  associated  with  correlated  hardware  errors  and  correlated  errors  due  to 
the  omission  of  linear  aggregate  deflection  terms.  These  latter  errors  are  not 
independent  from  one  traverse  to  an  adjacent  one. 
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5.  CONSIDERATION  OF  MOUNTAINOUS  TERRAIN.  The  use  of  advanced  or  Wiener-type 
methods  for  deflection  determination  in  semi-flat  terrain  requires  a  modifi¬ 
cation  in  mountainous  terrain  in  order  to  compute  signal  covariances.  As 
shown  by  Baussus  von  Luetzow  [1981b],  It  Is  possible  to  represent  deflections’ 
In  mountainous  terrain  in  the  form 

ft  -{!}♦{£} 

where  6£t  and  6nt  are  computable  topographic  "noise"  terms,  statistically  non¬ 
stationary  in  character.  It  is,  therefore  necessary  in  connection  with  the 
utilization  of  an  advanced  RGSS  to  employ  the  transformation 


q(2)|  =Ui-6Cti) 

ni(2))  (ni  -  6nt±) 


In  the  advanced  methods  discussed  in  sections  3  and  4.  This  results  in  the 
modification  of  measured  "message"  information. 

Following  these  data  modifications,  spatial  collocation  may  be  employed 
instead  of  planar  covariance  functions,  also  outlined  by  Baussus  von  Luetzow 
[1981b].  After  completion  of  the  signal  estimation  of  Its  correspond¬ 
ing  "message"  value  is 


Ce<2> 


*e  + 


<55 


te 


(21) 


Based  on  the  above,  the  optimal  estimation  of  deflections  of  the  vertical 
from  astrogeodetic  and  Inertial  data  is  not  possible  without  a  considerable 
computational  effort.  For  this  reason,  the  basic  method  presented  In  section  2 
appears  to  be  sufficient  If  data  from  one  or  two  repetitive  runs  are  employed. 
Regardless  of  these  repetitions,  the  deflections  accuracy  achieved  thereby  will 
be  somewhat  impaired. 
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6.  AREA  ADJUSTMENT  SOLUTIONS.  The  advanced  methods  presented  in  sections  4 
and  5  simultaneously  provide  for  area  adjustment  solutions.  To  assure  quasi- 
uniform  coverage  and  to  avoid  computational  complexity  and  the  Influence  of 
accelerometer  scale  factor  variations,  the  traverses  should  be  quasi -parallel 
and  system  calibration  should  be  attempted  prior  to  each  run  to  eliminate  cor¬ 
relations  between  systems  errors  referring  to  different  traverses.  Cross¬ 
traverse  coverage  would  be  beneficial  as  to  data  comparison  at  intersection 
points  and  for  additional  error  reduction.  Quasi-rectangular  coverage  with 
astrogeodetic  data  on  the  boundary  would  provide  for  economic  survey  extensions. 
For  practical  purposes,  approximately  50  adjacent  "message"  data  need  to  be 
used  for  a  specific  minimum  error  variance  determination.  Computer  program 
inputs  are  bias-corrected  "message"  deflections,  boundary  deflections,  their 
coordinates,  and  systems  error  variances  and  covariances  associated  with  signal 
covariances.  The  refined  advanced  method  requires  a  considerably  greater  pro¬ 
gramming  and  input  effort.  Because  of  the  reasonable  restriction  to  about  50 
"message"  data,  it  is  not  necessary  to  employ  mixed  data,  i.e.,  £-,  n-,  and 
gravity  anomaly  data.  Cross-covariances  involving  these  variables  are  insig¬ 
nificant  for  shorter  distances.  Further,  the  presence  of  gyro  heading  sensitiv¬ 
ity  errors,  reduced  or  unreduced,  does  not  warrant  undue  complexity.  As  shown 
in  section  5,  solutions  affected  by  formidable  mountainous  terrain  require 
special  consideration  or  will  be  impaired,  respectively. 

An  inverse-space  domain  smoothing  involving  the  use  of  orthogonal  functions 
and  deflection  and  gravity  anomaly  data  within  a  large,  flat,  rectangular  area 
has  been  developed  by  Bose  [1980].  In  the  context  of  this  method,  measurement 
noise  Is  assumed  to  be  a  zero  mean  uncorrelated  process.  Apart  from  this  approxi¬ 
mation,  the  "message"  data  has  to  be  generated  In  a  regular  pattern  which  Is  often 
not  possible.  In  order  to  achieve  an  effective  smoothing  under  avoidance  of  Gibbs 
type  fluctuations,  it  is  necessary  to  analytically  approximate  measurements  for 
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each  elementary  rectangle  for  the  computation  of  coefficients  a  In  the  case 

®n 

of  large  m  and  n. 

The  Wiener-type  solutions  under  utilization  of  approximate  signal  covari¬ 
ance  functions  and  non-stationary,  partially  correlated  system  errors  and  the 
Bose  solution  have,  therefore,  advantages  and  disadvantages.  The  former  ones 
are  certainly  more  versatile  and  in  principle  minimum  variance  estimates,  and 
the  refined  Wiener-type  solution  offers  an  additional  advantage. 


7.  CONCLUSION:  Advanced,  Wiener-type  filtering  methods,  including  refined 
methods  and  consideration  of  computable  topographic  "noise"  in  mountainous  areas, 
presented  in  this  paper,  will  provide  minimum  variance  estimates  of  deflections 
of  the  vertical  from  astrogeodetic  and  discrete,  inertially  derived  data  both 
for  single  and  multiple  traverses.  Additional  research,  to  include  numerical 
integration  of  a  system  of  differential  equations  and  Kalman  filter  corrections, 
is  necessary  to  facilitate  a  quantitative  comparison  between  the  method  of 
section  3  and  the  refined  method  of  section  4. 
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